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ABSTRACT 


This  report  presents  a  general  investigation  of  buffered  shock  t  bes. 

For  the  case  of  strong  shocks  the  number  of  variables  needed  for  the  analysis 
of  buffered  shock  tubes  is  reduced  to  three.  An  expression  is  obtained  from 
which  the  buffer  conditions  necessary  for  the  production  of  the  maximun  shock 
Mach  number  are  determined.  The  strong  shock  analysis  is  extended  to  include 
shock  tubes  with  many  buffers  and  buffered  shock  tubes  with  area  changes  at 
either  or  both  diaphragms.  The  concept  of  using  a  detonable  mixture  as  a 
buffer  gas  is  explored  and  its  limitations  are  indicated.  Also,  a  study  is  in¬ 
cluded  of  different  shock  tube  geometrical  configurations  to  indicate  means  of 
achieving  the  maximum  testing  time. 

The  strong  shock  assumption  implies  that  for  given  over-all  conditions 
the  final  shock  strength  depends  only  upon  the  upstream  density  ratio.  This 
simplification  permits  the  study  of  the  buffered  tube  throughout  its  entire  range 
of  operating  conditions.  A  number  of  general  conclusions  can  be  drawn  from 
this  study.  For  given  over-all  conditions  a  buffered  shock  tube  will  produce 
stronger  shocks  than  a  simple  shock  tube  provided  that  the  optimum  buffer  con¬ 
ditions  are  used.  The  performance  of  a  shock  tube  improves  as  the  number 
of  buffers  increases,  but  the  final  shock  strength  soon  approaches  an  asymptotic 
value.  The  performance  of  a  buffered  tube  is  improved  by  using  an  over-all 
area  contraction;  the  most  efficient  distribution  involves  equal  area  changes 
at  each  diaphragm.  The  gain  in  performance  obtained  with  a  detonation  buffer 
is  limited  when  the  final  shocks  are  very  strong.  For  fixed  driver  conditions 
and  final  shock  strength,  increasing  the  diameter  of  the  driven  section  does 
not  significantly  increase  the  testing  time. 
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1.  INTRODUCTION 


The  shock  tube  has  become  a  common  laboratory  tool  for  the  study  of 
aerodynamics,  chemical  kinetics,  combustion  processes,  and  high  tempera¬ 
ture  gas  physics.  The  literature  on  shock  tubes  is  very  extensive,  compilations 
of  published  works  can  be  found,  for  example,  in  Refs.  1  and  2.  In  recent  years 
the  shock  tube  has  been  extensively  used  as  a  driver  for  shock  tunnels  in  ■which 
Mach  numbers  and  stagnation  enthalpies  of  hypersonic  flight  can  be  simulated 
(Ref.  3  has  a  large  bibliography  on  the  subject).  Among  the  driver  gases  which 
have  been  used  for  these  devices  are  the  following;  cold  hydrogen,  cold  helium, 
hydrogen-oxygen  mixtures  with  and  without  diluent,  heated  hydrogen,  and  arc- 
heated  helium.  The  present  desire  to  simulate  flight  at  escape  velocity 
(36,000  ft/sec)  and  higher  has  stimulated  interest  in  extending  driver  cap¬ 
abilities.^’  ^ 

The  purpose  of  this  report  is  to  study  a  particular  method  for  increasing 

the  shock  Mach  number  in  a  shock  tube  with  any  given  driver.  This  method  is 

the  use  of  the  buffered  (or  multiple -diaphragm)  shock  tube.  The  use  of  buffer- 

16-18 

ing  has  been  considered  by  many  authors  during  the  past  ten  years.  ’ 

However,  because  of  the  complexity  of  the  problem,  the  previous  studies  have 
dealt  with  specific  examples  only,  and,  therefore,  do  not  reveal  in  a  general 
f  rm  the  maximum  gains  to  be  obtained.  The  present  report  utilizes  an 
i  sumption  cf  strong  shocks  to  simplify  the  analysis  of  optimum  performance 
With  this  simplification  the  optimum  buffer  conditions  are  predicted  over  the 
<  ntire  range  of  over-all  (i.e.  drive.*  to  driven)  conditions.  Analytical  expres 
smns  are  derived  for  the  maximum  final  shock  strength  obtainable  with  the 
optimum  buffer  conditions. 

Manuscript  released  by  the  author  Jan  iary.  1964  for  publication  as  an  ARL 
Technical  Documentary  Report. 
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It  may  be  noted  that  a  number  of  buffered  tubes  have  been  built, 
and,  in  general,  the  experimental  results  support  the  theory.  However,  the 
actual  performance  of  these  tubes  is  below  that  predicted  by  ideal  theory  due 
to  such  effects  as  shock  wave  attenuation,  wave  reflections,  contact  zone 
mixing,  and  real  gas  effects. 

The  convenience  of  using  strong  shock  assumptions  to  simplify  the  pre¬ 
sent  analysis  of  buffered  shock  tubes  was  pointed  out  in  the  work  of  Hall  and 
19 

Russo.  They  present  a  useful  simplification  to  the  simple  shock  tube  equation 
for  shock  Mach  numbers  exceeding  about  three.  For  a  given  driver  gas  specific 
heat  ratio,  shock  tube  performance  can  be  explicitly  expressed  for  all  initial 
conditions  (including  area  change)  in  terms  of  two  variables.  They  note  that 
the  performance  of  buffered  tubes  (with  strong  shocks)  can  be  expressed  in 
terms  of  corresponding  variables  with  only  the  upstream  dc..  _  ty  ratio  as  an 
additional  variable.  They  further  state  that  previous  results  of  Russo  and 
Hertzberg*^  showing  the  effects  of  varying  buffer  gas  pressure  and  molecular 
weight  can  be  correlated  in  terms  of  this  density  ratio. 

In  Section  2.  1  a  survey  of  previous  studies  of  buffered  shock  tubes  is 
presented.  The  theoretical  studies  all  bear  out  the  fact  that  the  performance 
of  any  shock  tube  can  be  improved  by  the  inclusion  of  a  buffer  section.  Cal¬ 
culations  of  the  performance  of  a  buffered  shock  tube  for  any  given  initial 
conditions  in  the  driver  and  driven  sections  are  quite  lengthy  because  of  the 
number  of  variables  involved:  buffer  gas  specific  heat  ratio,  molecular 
weight,  temperature,  and  pressure.  For  a  given  buffer  gas  at  a  fixed  initial 
temperature  the  final  shock  Mach  number  varies  only  with  the  buffer  pressure. 
There  exists  a  certain  buffer  gas  pressure  (which  shall  be  called  the  optimum 
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pressure)  for  which  the  final  shock  strength  is  a  maximum.  The  determina¬ 
tion  of  this  optimum  pressure  is  complicated  by  the  fact  that  the  calculation 
ol  the  driven  section  shock  Mach  number  involves  two  iterations  for  each 
set  of  initial  conditions. 

In  Section  2.  2  the  equations  are  simplified  with  the  assumption  of 
strong  shock  waves  so  that  analytical  expressions  for  the  optimum  buffer  con¬ 
ditions  are  obtained.  V/hen  the  optimum  conditions  have  been  determined, 
analytical  expressions  are  obtained  for  the  maximum  final  shock  strength  for 
any  given  over-all  conditions.  The  strong  shock  analysis  is  extended  with 
similar  results  to  the  study  of  shock  tubes  with  many  buffers  and  to  buffered 

tubes  with  area  changes  at  either  or  both  diaphragms. 

12 

The  use  of  a  detonable  mixture  as  a  buffer  gas  is  considered  in 
Section  2.  3.  This  device  has  most  of  the  characteristics  of  a  buffered  shock 
tube  with  the  added  advantage  that  chemical  energy  is  released  to  the  flow  in 
the  buffer  section.  When  the  strong  shock  assumption  is  applied  to  the  detona¬ 
tion  buffer  equations  (assuming  ideal  gases),  the  optimum  performance  con¬ 
ditions  can  be  determined.  If  the  tube  is  operated  so  that  the  detonation  in 
the  buffer  is  close  to  a  Chapman- Jouguet  detonation,  a  substantial  increase  in 
final  shock  Mach  number  is  possible.  However,  if  the  detonation  is  strongly 
ov^r-driven,  the  detonation  buffer  offers  only  moderate  gains  over  the  inert- 
gas  buffer. 

In  Section  3  studies  are  made  to  determine  the  geometrical  configuration 
which  yields  the  greatest  testing  time  for  a  given  final  shock  strength  and 
specified  driver  conditions.  It  has  been  pointed  out  that  there  is  a  severe  de¬ 
crease  in  flow  duration  for  shock  tubes  operating  at  low  pressures  due  to  the 
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The  testing  time 


effects  of  the  boundary  layer  behind  the  shock  waves, 
as  limited  by  boundary  layer  growth  has  been  shown  to  be  directly  proportional 
to  the  pressure  and  to  the  area  of  the  driven  section.  This  fact  has  stimulated 
interest  in  shock  tubes  with  larger  driven  sections  than  driver  sections.  For 
most  cases  of  interest  the  benefit  of  the  larger  diameter  is  found  to  be  offset 
by  the  need  for  a  lower  initial  pressure  to  maintain  the  same  shock  Mach 
number. 

2.  OPTIMUM  PERFORMANCE  OF  BUF  FERED  SHOCK  TUBES 
2.  1  REVIEW  OF  PREVIOUS  ANALYSES 

The  use  of  the  double-diaphragm  shock  tube  (i.  e.  the  buffered  shock 

1  o-18 

tube)  has  been  considered  by  a  number  of  authors.  ’  In  the  following 

analyses  it  will  be  assumed  that  the  flow  is  one-dimensional  and  inviscid  with 
no  heat  conduction  and  no  mixing  at  the  interfaces.  For  the  most  pari  the  gases 
will  be  assumed  to  be  ideal  gases  with  constant  specifi  :  heats. 

The  basic  arrangement  is  illustrated  in  Figure  1  for  the  "unsteady  - 

expansion"  type  of  buffered  shock  tube.  The  initial  flow  phenomena  are 

represented  on  a  distance  vs.  time  (x-t)  plot  which  is  commonly  known  as  a 

wave  diagram.  The  high  pressure  in  the  driver  (region  -!)  causes  the  initial 

— ^ 

diaphragm  to  rupture.  A  shock  wave  5,  propagates  into  the  buffer  gas, 
while  a  rarefaction  (or  expansion)  wave  ^  propagates  into  the  high  pressure 
driver  gas.  The  driver  gas  and  the  buffer  gas  remain  separated  by  a  contact 
surface  across  which  the  pressure  and  velocity  are  equal,  but  the  density, 
temperature,  and  entropy  are,  in  general,  different. 

Assuming  that  the  second  diaphragm  bursts  instantaneously  upon  the 
impact  of  the  initial  shock,  the  wave  diagram  is  represent*  d  by  Figure  1.  The 
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flew  from  region  2  to  region  8  resembles  that  of  a  simple  shock  tube  except 
that  the  "driver"  has  an  initial  velocity  U%  and  has  been  heated  by  the  shock 
to  a  temperature  .  However,  the  pressure  is  less  than  the  initial 

driver  pressure  'P'+-  Therefore,  the  double-diaphragm  shock  tube  can  be 
considered  as  a  device  which  effectively  sacrifices  pressure  ratio  in  favor 
of  an  increased  speed  -  of- sound  (or  temperature)  ratio. 

A  second  possibility,  termed  the  reflected- shock  type  of  buffered  shock 

1  9 

tube,  is  illustrated  in  Figure  2.  ’  With  this  configuration  the  initial  shock  is 
allowed  to  undergo  complete  reflection  from  the  second  diaphragm  leaving  the 
gas  in  region  5  at  rest.  After  a  predetermined  delay  Y  the  second  diaphragm 
is  ruptured  mechanically;  the  ensuing  flow  is  shown  in  Figure  2.  While  this 
configuration  offers  a  slight  improvement  over  the  unsteady- expansion  type, 
the  mechanical  difficulties  involved  in  breaking  the  second  diaphragm  will 
preclude  its  further  consideration.  The  reflected- shock  type  buffered  tube 
should  net  be  confused  with  the  shock  tube  with  an  area  contraction  at  the 
second  diaphragm.  This  tube  also  produces  a  reflected  shock,  but  the  ensuing 
flow  is  fundamentally  different  from  that  of  the  reflected- shock  type  of  tube. 

Consider  some  of  the  characteristics  of  the  unsteady-expansion  buffered 
shock  tube.  Before  the  shock  strikes  the  second  diaphragm  the  flow  is  identical 
to  that  which  is  found  in  a  simple  shock  tube.  Therefore,  the  well-known 
simple  shock  tube  equation  is  valid  (see  Appendix  I). 


2?-,  <  -(r,-Q 

Z+l 


7 K  -  ') 

T‘+'  \,MS, 


V' 


(1) 
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where 


initial  shock  Mach  number 


'  Vs,/*! 

L/  speed  of  the  shock  relative  to  the  tube 

(L  speed  of  sound 

p^(  =  ^>. initial  diaphragm  pressure  ratio 
A+(=  initial  diaphragm  speed-of- sound  ratio 

/  r  ratio  of  specific  heats 

Referring  to  Figure  1  and  applying  the  usual  equations  of  one-dimensional 
unsteady  flow,  the  basic  equation  for  the  constant-area  buffered  shock  tube  is 
obtained  (see  Appendix  I). 


where  Ms  final  shock  Mach  number 

n4|  =  p  Jp  density  ratio 

The  final  shock  strength  is  a  function  of  the  following  eight 

8 

variables:  ^  ,  p  ,  p  •  M,  '  anc*  •  However,  these  ar 

not  all  independent,  and  one  can  be  eliminated  by  the  use  of  the  simple  shock 

tube  equation  (Eq.  (l)  ).  As  in  the  case  of  the  simple  snock  tube,  M,  is 

increased  if  p  is  made  large  and  P.„  is  made  small  (i.e.  making  /[ 

48  °  +8 

large).  Also,  as  with  the  simple  shock  tube,  the  best  performance  is  obtained 
when  7^  is  small  which  suggests  a  diatomic  driver  gas.  The  effect  of 


f> 


varying  y  is  found  to  be  negligible;  a  monatomic  driven  gas  is  found  to  be 
only  slightly  better  than  a  diatomic  gas. 

When  conditions  in  the  driver  and  the  driven  sections  (4  and  8)  have  been 
specified,  M  depends  upon  the  three  variables  y  ,  P  ,  and  F  .  Russo 

S  ^  *  Tl  tI 

and  Hertzberg*^  have  shown  that  the  performance  of  a  buffered  shock  tube 
improves  with  increasing  buffer  specific  heat  ratio  jr  .  Also,  the  use  of 
monatomic  buffer  gases  will  minimize  real  gas  effects  which  tend  to  reduce 
performance . 

9 

Henshall  presents  results  for  a  buffered  tube  in  which  the  same  diatomic 
gas  at  the  same  initial  temperature  is  used  throughout.  In  this  case  yf  -  y^  r 
7a  -  7/5,  P^  =  and  yU,=  where  yU  is  the 

molecular  weight.  The  final  shock  Mach  number  depends  only  upon  the  over -all 
pressure  ratio  and  the  upstream  pressure  ratio.  Figure  3  shows  as  a 

function  of  the  downstream  pressure  ratio  Pjg  with  P ^  as  a  parameter.  It 
will  be  noted  that  the  maximum  value  of  M.  occurs  when  P|fl  is  approximately 


A  similar  result  can  be  seen  from  Figure  4  which  shows  the  per¬ 


formance  of  a  tube  in  which  the  same  monatomic  gas  at  the  same  initial  temper¬ 
ature  is  used  throughout.  Of  course,  in  a  more  general  case,  the  geometric 
mean  may  not  be  the  optimum  buffer  pressure.  The  shock  Mach  numbers 
obtained  from  simple  shock  tubes  with  the  same  over-all  conditions  are  in¬ 


cluded  in  Figures  3  and  4  for  comparison. 

Hydrogen  is  a  widely  used  driver  gas  since  it  has  tne  smallest  molecular 

weight,  anc  therefore,  the  smallest  density  at  a  given  temperature  and  pressure. 

Figure  5  shows  the  performance  of  a  tube  with  hydrogen  driver,  monatomic 

4 

buffer,  and  -lir  as  the  driven  gas  at  an  over-all  pressure  ratio  of  1  0  .  The 


7 


temperature  is  assumed  constant  throughout  so  that 

rf8  =  P48A,  and  O) 

where  /M4g=^4^^*  the  ratio  of  molecular  weight;-.  Therefore,  is  a 

fi 

function  only  of  U,  and  P 

A  study  of  Figure  5  reveals  the  following  .nteresting  feature,  while  the 

optimum  downstream  pressure  ratio  P  vines  with  the  buffer  gas  molecular 

weight,  the  maximum  value  of  M,  is  relatively  unaffected.  This  seems  to 

8 

indicate  that  certain  simplifications  of  the  basic  equations  are  possible  which 
reduce  the  number  of  independent  variables.  It  will  be  shown  presently  that 
this  is  indeed  the  case  and  that  the  simplified  buffer  equation  may  be  used  to 
determine  the  optimum  buffer  conditions. 

2.  2  INERT-GAS-BUFFERED  SHOCK  TUBES 


2.2.1  Simple  Constant-Area  Buffered  Tubes 

1  9 

Recently  Hall  and  Russo  have  shown  that  the  analysis  of  a  simple 
shock  tube  with  an  area  change  near  the  diaphragm  can  be  greatly  simplified 
by  assuming  that  the  shock  Mach  number  is  greater  than  about  three.  This 


theory  is  reviewed  in  Appendix  II  where  it  is  shown  that  for  given  7  and  7 

M  ' 

the  basic  equation  involves  only  two  variables:  j?'  =  y  and  P^  -  ^ 

where  is  the  equivalence  factor  of  Rosier,  Lin,  and  Kantrowitz.^  A  plot 

of  the  basic  relation  is  given  in  Figure  6  for  various  values  of  7^  the 
variation  in  performance  with  is  found  to  be  negligible.  The  plot  can  be 

interpreted  as  representing  the  performance  of  a  simple  shock  tube  merely  by 


setting  2*  equal  to  one. 

Now  apply  the  strong  shock  assumptions  to  the  constant-area  buffered 


8 


tube.  Assuming  that  both  and  M  are  greater  than  about  three,  Eq.  (2) 

5 1  ’a 

can  be  simplified  to  the  relation 


This  equation  can  be  further  simplified  by  the  use  of  the  simple  shock  tube 
equation,  which  in  this  case  becomes 


Equation  (5)  is  easily  solved  for  P^  yielding 
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where 


Y.  +  1 

c.  -  6 


Y*l 


For  given  specific  heat  ratios  and  fixed  over-ali  density  ratio  P  a  ,  ZB 

0  ifl  5 

deoends  utx>n  K  and  2  .  It  is  possible  to  maximize  Z  with  respect  to 

*,  '  8 

these  two  variables  and  obtain  the  conditions  for  optimum  performance.  Ihe 
details  of  this  calculation  are  presented  in  Appendix  III.  The  important  point 
is  that,  once  the  maximizing  is  done,  the  final  shock  strength  £  is  a  function 


only  of  the  over-all  density  ratio  f 


•*8 


Upon  now  making  the  additional  assumption  that  NJ  is  large  compared 


7,-1 


,  Eq.  {7}  reduces  to  a  particularly  simple  form. 


.  J 

i  ~  r  ?-  J  D  - 

8  1  * 


8£,P  '/a  Z 

•  +£  8 


(fr  ") 


(8) 


where 


and 


D,  =  I  + 


/Vi_f 
\  2/,/ 


p-'f , 

u/  r,*t 

Equation  {8}  can  be  -written  symbolically  as 

zj  = F.  fr  *  .  'i .  ■  r«) 


(9) 


When  the  over -all  density  ratio  and  the  specific  heat  ratios  of  all  gases  are 
specified,  the  final  shock  strength  depends  only  upon  Z,  •  Alternatively, 

Eq.  (6)  may  be  used  to  eliminate  the  explicit  dependence  of  Zg  upon  Z ; 

Then,  for  specified  over-ali  density  ratio  and  specific  heat  ratios,  Z  would 
depend  only  upon  the  upstream  density  ratio  f1  .  This  result  was  previously 
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(10) 


noted  in  Ref.  19.  The  value  of  2  needed  to  yield  a  maximum  2  should 

*  s 

be  a  solution  to  the  equation 


32. 


-0 


The  details  of  this  calculation  are  shown  in  Appendix  III  where  it  is  found 


that 


and 


P48  = 


H e# 


B  E,  D,  C,(l+x)  Z, 


8 


■(W 


(ID 


(12) 


where 


X,  = 


(13) 


Equations  (11)  and  (12)  may  be  used  to  determine  2  —maximum  as  a 

O 

function  of  by  means  of  the  parametric  dependence  of  these  two  quantities 
upon  2  .  Thus,  for  any  given  P  ,  2  (optimum)  can  be  calculated  from 

•  TO  • 

(11),  the  required  upstream  density  ratio  P  can  be  found  from  (6),  and  the 
maximum  final  shock  strength  can  be  obtained  from  (12).  Figure  7  is  a  plot 
of  Eq.  (11)  for  the  special  case  in  which  y  -  -  1  IS  and  Y ,  -  SIS,  while 

Fig.  8  of  2-  vs.  r.  Q  is  obtained  from  Eqs.  (11)  and  (12).  Also,  on  Fig.  8, 

0*0 

for  comparison,  is  a  plot  of  2,  vs.  P  for  a  simple  shock  tube  w'  h  Y^  - 


As  an  example  of  the  calcuiational  procedure,  consider  a  tube  with  a 
hydrogen  driver,  a  monatomic  buffer,  and  driven  air  all  initially  at  the  same 
temperature.  Then 


P4»  =  P« 

—  4 

In  order  to  compare  this  with  the  results  of  Fig.  5  let  =  10  ,  P. .  =  696. 


Then  from  Fig.  7  =  .189,  from  Fig.  8  2F  =  .117,  and  from  Fig.  6 

P  =  130.  Now  since  R.  =  10  ,  M*  =11.7  which  is  in  good  agreement 

♦t  *  g 

with  Fig.  5.  If  we  assume  that  the  buffer  gas  is  argon  (  =  39.9)  then 

P+l  ~  *4i/^4|  =  2580  or  Pi9  =  3.88  which  is  again  in  excellent  agreement 

with  Fig.  5. 

Some  idea  of  the  validity  of  the  assumptions  can  be  obtained  from  Fig.  9 
which  shows  a  comparison  of  the  exact  solution  (Eq.  (2)  ),  approximation  1 
(Eq.  {"?)  ),  and  approximation  2  (Eq.  (8)  )  for  two  different  monatomic  buffer 
gases  -  4  and  y/U.t  80).  The  plots  for  other  values  of  the  molecular  weight 

are  similar.  Approximation  2  becomes  nolicably  worse  as  Plg  becomes  large; 
this  is  to  be  expected  since  Ms  will  be  small  for  large  Pf^  .  Since  the  maxi¬ 
mum  f°r  the  ll  =80  gas  occurs  at  small  p  ,  approximation  2 

reliably  predicts  the  maximum  while  the  approximation  for  the  yU  f  =4  case 
is  less  accurate.  This  and  other  calculations  indicate  that  approximation  2 

alwjvs  yields  values  cf  M.  lower  than  the  exact  solution,  so  that  it  can  be  con- 

7  s8 

sidered  as  a  lower  bound.  The  curves  of  M5  vs.  P.8  for  approximation  2  are 

8 

all  txactly  similar  (no  matter  what  the  buffer  gas  molecular  weight),  while  the 
curves  for  approximation  1  show  a  slight  dependence  of  2  max.  upon 

O 

molecular  weight  (just  as  the  exact  curves  do). 

P'igure  10  presents  another  comparison  of  the  approximations,  this  time 


1  2 


for  the  case  where  -  y  =  7/5,  -  /l%  .  /,  =  7^  =  Tg  and 

P  =  lO^.  In  this  case  approximation  2  gives  a  value  of  Ms  which  is  in 

error  by  about  15%  near  the  optimum.  This  error  is  due  to  the  low  values 

of  M  (around  3)  which  occur  near  the  optimum  conditions.  When  the  pro- 

duction  of  strong  shocks  is  a  primary  concern,  however,  the  values  of  M«. 

*'1 

which  are  needed  are  usually  such  as  to  make  approximation  2  reasonably  good. 

In  Appendix  III  the  problem  of  determining  the  optimum  conditions  from 
approximation  1  is  given  detailed  consideration.  No  analytical  expression  was 
found  from  which  the  optimum  conditions  for  any  could  be  determined. 

Solutions  to  the  problem  were  obtained  for  fixed  values  of  P^.  and  K 


where 


ami 


=  '  41  _  >^4I 


41 


T 


♦i 


T  -  Jj— 

•w  ~  r, 

but  the  calculations  are  quite  lengthy.  However,  the  following  important  fact 
can  be  noted:  for  fixed  P  the  value  of  —maximum  increases  as  the 

value  of  increases.  This  result  is  illustrated  in  Fig.  9  which  shows 

that  £8— maximum  is  greater  for  yU/  =  4  (  K*  =  .5)  than  for  yU  t  -  80 
(  fc  =  .  025).  It  is  interesting  to  note  that  as  K  increases,  the  value  of 

at  —maximum  decreases,  and,  therefore,  the  strong  shock  assump¬ 

tion  is  becoming  increasingly  invalid. 

The  difficulties  involved  in  determining  the  optimum  conditions  with  the 
use  of  approximation  1  are  comparable  to  those  involved  with  the  use  of  the 
exact  equation.  The  calculations  of  the  optimum  conditions  with  approximation 
2,  on  the  other  hand,  are  relatively  straightforward.  Even  in  cases  where  the 
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second  approximation  is  not  very  accurate,  approximate  values  for  the 
optimum  conditions  can  be  obtained  for  use  as  initial  guesses  for  exact 
calculations. 

2.2.2  The  Many -Buffered  Tube 

A  shock  tube  with  more  than  one  buffer  section  should  produce  a  stronger 

9 

final  shock  than  a  simple  buffered  tube.  Henshall  considers  the  special  case 
of  all  diatomic  gases  with  the  same  molecular  weight  and  initial  temperatures. 

He  assumes  that  the  optimum  pressure  ratio  distribution  will  be  equal  pressure 
ratios  across  each  diaphragm.  Under  these  assumptions  the  final  shock  strength 
increases  with  the  number  of  diaphragms  but  the  relative  gain  decreases  as  the 
number  of  diaphragms  increases.  As  the  number  of  diaphragms  becomes  very 
large,  the  final  shock  Mach  number  approaches  an  asymptotic  value. 

The  calculation  of  the  optimum  conditions  for  a  many -buffered  tube  is 
very  difficult  because  the  final  shock  strength  depends  upon  the  density  and  the 
pressure  in  every  section  of  the  tube.  It  is  therefore  very  desirable  that  a 
method  be  developed  to  determine  quickly  these  optimum  conditions.  The 
basic  equation  for  the  many-buffered  tube  is  developed  in  Appendix  I  and  is 
based  upon  an  idealized  wave  diagram  such  as  that  shown  in  Fig.  1 1 . 

The  wave  diagram  in  Fig.  11  indicates  the  numbering  system  that  will  be 
used  when  discussing  the  many-buffered  shock  tube.  The  system  used  for 
numbering  the  flow  regions  in  simple  shock  tubes  (driver-4,  driven- 1)  is  almost 
universally  accepted.  The  numbering  system  for  the  buffered  shock  tube  which 
is  most  widely  used  is  that  which  labels  the  driven  section  8  and  the  region  be¬ 
hind  the  reflected  shock  (if  any)  5.  Figure  11  represents  the  extension  of  this 
system  to  a  double-buffered  tube.  Regions  5  and  9  are  omitted  because  there 


J  4 


are  no  reflected  shocks. 

For  simplicity,  begin  by  considering  a  double-buffered  shock  tube.  The 

basic  equations  are  (A- 23),  the  simple  buffer  Eq.  (2),  and  the  simple  shock 

t  a  a 

tube  Eq.  (1).  Assuming  that  M  ^  I  »  and  M_  >}  J  >  an  initial 

S'  s8  S»2 

approximation  is  obtained,  but  the  equations  are  still  so  complex  that  no 

attempt  will  be  made  to  maximize  them.  When  the  rather  restrictive  assump- 
z  2  2  2 

tions  that  M  - —  and  M  ^  - — -  are  made,  it  is  found  that  Eq.  (1)  reduces 

i  7i~  I 

to  Eq.  (5),  Eq.  (2)  reduces  to  Eq.  (8),  and  Eq.  (A-23)  reduces  to 


Now  Eq.  (8)  can  be  solved  for  p  yielding 

»  o 


1  '/* 
+8 


So  that  Eq.  (14)  becomes 


BE, 


(14) 


(15) 


15 


Equation  (16)  can  be  represented  symbolically  as 

^  12  ~  Fz  >  Zf  >  t  Kz  >  Zj  >  ^g,  r4  #a J  (17) 

So,  for  given  specific  heat  ratios  and  over-all  density  ratio,  the  conditions  for 
maximum  zL/z  are 


3?, 


=  0 


and 


(18) 


The  calculations  are  quite  lengthly,  but  fortunately  the  results  can  be 
expressed  as  equations  involving  only  the  single  independent  variable  . 
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Figure  12  shows  ?|2  vs.  P4>ia  calculated  from  Eqs.  (19)  and  (20),  for 

the  special  case  in  which  /(  =  7^  =  =  5/3.  This  result  is  shown 

compared  with  the  results  from  the  simple  shock  tube  equation  and  the  simple 

buffered  tube  equation.  The  double-buffered  tube  curve  lies  above  the  other 

(a.A  \ 

curves,  the  relative  improvement  increasing  near  the  two  limits  P^  (2 — >0\^ - 

and  P  ->oo (p  — Voo)  .  It  is  interesting  to  note  that  in  one  region  (Pm  11. 5, 

2  ~  .36)  no  improvement  is  obtained  by  using  buffers.  This  point  was  checked 
by  several  exact  calculations  arid  found  to  be  quite  accurate,  so  that  in  this 
range  a  simple  shock  tube  is  just  as  efficient  as  a  buffered  one.  It  will  also  be 
noticed  that  the  relative  gain  of  the  double-buffered  tube  over  the  simple 
buffered  tube  is  always  less  than  the  relative  gain  of  the  simple  buffered  tube 
over  the  simple  shock  tube.  Similar  results  are  indicated  in  Figs.  13  and  14 
for  different  gas  combinations. 

Evaluating  the  triple -buffered-tube  equations  (assuming  strong  shocks), 
it  is  found  that 


so  that  the  following  conditions  must  be  satisfied  to  determine  the  maximum 
final  shock  strength: 
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15.  =o  ,  15 


-  0 


The  final  results  are  quite  similar  to  Eqs.  (19),  (20)  and  (21).  In  fact,  the  resulting 
equations  indicate  a  general  trend,  so  that  it  is  possible  immediately  to  general¬ 
ize  to  a  many-buffe red  tube  of  any  number  of  sections.  For  example,  if  the 
driven  section  is  called  number  (71+4)  (based  on  the  general  numbering  system 
indicated  in  Fig.  11),  then 


//-#-  x 


v^;  c„  e»*o* 


2  =  D 


r „ 

V'  'A 


.  +  X, 


where 


rw*-i 


= 


r, 


D,  =  /  + 


r,-/VA 


(W 


V ' 
+ ' 


Cn  = 


y*f4-  +  1 


These  equations  are  valid  only  under  the  assumption  that  all  the  inter¬ 
mediate  shock  Mach  numbers  are  greater  thin  about  seven.  Of  course,  as  the 


number  of  sections  becomes  large,  these  assumptions  will  generally  be  violated. 
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The  results  obtained  from  (24)  and  (25)  are  then  only  qualitative  indications  of 
performance,  but  as  such  they  indicate  an  important  trend.  The  final  shock 
strength  rapidly  approaches  some  asymptotic  value  as  the  number  of  dia¬ 
phragms  increases,  and  this  value  depends  only  upon  the  over-all  density 
ratio.  In  actual  practice  the  performance  of  a  many-buffered  tube  will  depend 
critically  upon  the  lengths  chosen  for  each  section.  Wave  interactions  will 
generate  new  waves  which  travel  in  the  same  direction  as  the  final  shock  and 
tend  to  limit  testing  time. 

2.  2.  3  Effects  of  Area  Contractions 

A  method  for  improving  the  performance  of  a  simple  shock  tube  is  to 
employ  a  monotonic  area  contraction  near  the  diaphragm  (see  Appendix  II). 
The  performance  of  such  a  shock  tube  is  the  same  as  that  of  a  simple  shock 


tube  with  a  pressure  ratio  of  and  a  speed  of  sound  ratio  of  y.  *  A4J 

where  ^  is  the  equivalence  factor.  Obviously,  the  performance  of  a  buffered 
shock  tube  could  be  improved  by  the  use  of  an  area  contraction,  but  here  a 
contraction  may  be  put  at  the  first  or  second  diaphragm  or  both. 

Consider  a  contraction  at  the  first  diaphragm  (see  Fig.  15).  The  initial 
flow  from  region  4  to  region  1  is  just  the  same  as  that  in  a  simple  shock  tube 
with  an  area  change  at  the  diaphragm  (Eq.  (A-34)  ),  while  the  flow  from  region 
2  to  region  8  is  the  same  as  that  in  a  constant-area  buffered  tube.  Therefore, 
the  basic  equation  for  this  configuration  is 
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The  analysis  of  the  buffered  tube  with  the  area  contraction  at  the  second 
diaphragm  is  a  much  more  difficult  problem.  The  general  theory  is  developed 
in  Appendix  IV  where  it  is  shown  that  there  are  three  possible  wave  diagrams 
when  M  v>/  (  see  Fig.  16).  The  case  of  most  practical  interest  is  the  one  for 

which  the  area  contraction  will  be  great  enough  to  produce  a  reflected  shock 
(Fig.  16-c).  Under  these  conditions  the  flow  is  adequately  represented  by 
Eas.  (A-72)  and  (A-76).  Unfortunately,  (A-76)  is  a  fourth  order  algebraic 
equation  for  the  reflected  shock  Mach  number,  and  it  is  usually  solved 

by  iteration. 

If  the  shock  tube  has  area  changes  at  both  diaphragms,  the  problem  can 

be  analyzed  by  a  combination  of  the  above  two  methods.  Figure  17  shows  the 

performance  of  a  He~a^r  shock  tube  operating  at  constant  initial  temper- 

4 

ature  and  with  an  over-all  pressure  ratio  of  1  0  .  When  the  over-all  contrac¬ 
tion  ratio  is  nine,  Fig.  17  shows  the  effect  of  different  distributions  of 
contraction  ratios.  The  best  performance  is  obtained  with  equal  area  ratios 
at  each  diaphragm.  If,  however,  it  is  desired  to  have  all  of  the  contraction 
at  one  diaphragm,  the  better  performance  is  obtained  with  the  area  change  at 
the  second  diaphragm. 
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When  the  area  change  is  all  at  the  second  diaphragm,  an  analytical  solu¬ 
tion  to  the  optimization  problem  can  be  ODtained  only  if  the  contraction  is 
infinite.  Fortunately,  for  area  ratios  of  ten  or  greater,  the  area  change  may 
be  considered  essentially  infinite.  In  this  case 


fa-OM,1  t  Z 

•f 


(31) 
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and  the  basic  Eq. 
assumption  that 

9  7,-1 

where 


(A-72)  is  greatly  simplified  (see  (A-78)  ).  When  the 


and 


/ 


(32) 


Qrm  is  the  equivalence  factor  for  an  infinite  area  contraction  at  the  second 
diaphragm.  Obviously  can  be  maximized  with  respect  to  in  exactly 
the  same  way  as  for  a  constant-area  buffered  tube.  The  results  are  as  follows: 


(34) 


In  Fig.  18  a  comparison  is  shown  between  the  exact  Eq.  (A-78)  and  the  strong 
shock  Eq.  (32)  for  the  case  of  cold  hydrogen  driver,  monatomic  buffer,  and 
driven  air.  The  agreement  is  excellent  at  large  values  of  buffer  gas  molecular 
weight  and  reasonably  good  even  for  a  helium  buffer.  Figure  18  is  entirely 
analogous  to  Fig.  9  except  for  the  area  change,  and  the  discussion  of  Fig.  9 
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applies  equally  well  here. 

Figure  19  can  be  constructed  from  Eqs.  (33)  and  (34).  The  graph  again 
points  up  the  fact  that  some  improvement  may  be  obtained  by  using  gases  with 
small  values  of  ,  but  the  performance  is  quite  insensitive  to  changes  in  . 

2.  3  DETONATION  BUFFERS 


One  method  of  improving  shock  tube  performance  is  the  use  of  some  com- 

1  (y  1  i  1 7  24 

bustion  process  in  the  driver  to  raise  the  pressure  and  temperature.  ’  ’  »  -  » 

The  most  widely  used  of  these  techniques  is  constant- volume  combustion  in  which 
the  gas  is  ignited  and  the  diaphragm  ruptures  just  as  the  peak  pressure  is 
reached.  Performance  calculations  for  this  type  of  shock  tube  can  be  made  by 
using  the  simple  shock  tube  equation  and  effective  values  of  ,  7^  ,  and  a 
corresponding  to  conditions  after  combustion.  Therefore,  the  previous  methods 
of  analysis  are  applicable  to  any  shock  tube  with  a  constant-volume  combustion 


driver. 


Two  other  combustion  processes  which  have  been  used  for  driver  gas 
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heating  are  constant-pressure  combustion  and  detonation.  Both  of  these 
methods  suffer  because  the  initial  shocks  are  closely  followed  by  rarefaction 
waves  which  limit  steady  flow  and  cause  shock  wave  attenuation.  Also,  there  is 
a  lack  of  uniformity  of  diaphragm  ruptures  with  constant -pressure  combustion. 

The  fundamental  equations  of  detonation  wave  theory  are  reviewed  in 
Appendix  V  where  it  is  shown  that  a  pure  detonation  is  always  followed  by  an 
expansion  wave  that  lowers  the  pressure  and  velocity.  If  a  piston  were  main¬ 
tained  in  motion  behind  the  detonation  with  a  velocity  equal  to  that  of  the  flow 
immediately  behind  the  detonation  front,  no  expansion  could  exist.  If  the  piston 
where  maintained  at  a  higher  velocity,  the  detonation  wave  would  be  propagated 
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at  a  faster  speed,  that  is,  it  would  be  overdriven. 

Consider  a  constant-area  buffered  tube  with  a  detonable  mixture  in  the 
buffer.  Such  a  device  will  be  called  a  detonation  buffer.  If  the  conditions  in 
the  driver  section  are  such  that  and  are  equal  to  or  greater  than 

the  Chapman-Jouguet  values,  the  detonation  in  the  buffer  will  be  followed 
by  a  region  of  quasi- steady  flow.  The  wave  diagram  will  then  be  similar  to 
that  shown  in  Fig.  1,  but  the  calculations  will  not  be  so  straightforward.  In 
particular,  the  flow  through  the  unsteady  expansion  between  regions  2  and  6 
must  be  calculated  by  integrating  the  characteristic  equation 

*«=  -f  V 


(35) 


along  an  isentrope.  ^ 

Some  exact  calculations  were  performed  for  a  shock  tube  with  a  hydrogen 
driver  gas,  air  driven  gas,  and  a  stoichiometric  hydrogen- oxygen  buffer  gas 
(with  an  without  helium  diluent).  Equation  (35)  was  integrated  with  the  help  of 
Mollier  charts  constructed  from  data  obtained  from  Dr.  Sanford  Gordon  at  JLewis 
Research  Center,  NASA.  Conditions  behind  the  over-driven  detonation  were 
found  by  using  the  equilibrium  normal  shock  IBM  704  program  developed  by  Dr. 
Russell  E.  Duff  of  Los  Alamos  Scientific  Laboratories.  Throughout  the  calcu¬ 
lations  the  hydrogen  driver  was  considered  as  an  ideal  gas,  and  the  air  was 
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treated. as  being  in  equilibrium.  '  The  exact  calculations  are  very  time 
consuming,  and,  consequently,  only  a  few  of  them  were  made.  Some  typical 
results  are  shown  in  Fig.  20  as  a  comparison  to  the  approximate  theory  that 
is  developed  below. 

The  addition  of  helium  diluent  to  stoichiometric,  hydrogen-oxygen  has 

28 

some  interesting  effects.  As  helium  is  added,  the  detonation  velocity 
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increases  and  the  molecular  weight  of  the  burnt  gas  decreases.  The  pressure 
and  temperature  behind  the  detonation  drop,  but  the  speed  of  sound  increases. 

As  a  net  result,  the  final  shock  strength  does  not  vary  much  with  the  addition 
of  helium.  Of  course,  the  detonation  buffer  does  offer  an  improvement  in 
performance  over  the  inert  gas  buffer,  but  the  gain  is  not  outstanding. 

In  order  to  compare  the  detonation  buffer  with  the  constant-area  inert- 
gas-buffered  tubes  over  the  full  range  of  operating  conditions,  it  is  desirable 
to  maximize  the  final  shock  strength  and  present  a  general  vs.  7  curve 
for  the  detonation  buffer.  To  facilitate  such  an  analysis  a  number  of  assumptions 
will  have  to  be  made.  Suppose  that  the  detonable  mixture  may  be  treated  as  an 
ideal  gas  before  detonation  and  as  a  different  ideal  gas  after  detonation.  Further 
assume  that  the  unsteady-expansion  flow  can  be  calculated  in  terms  of  an 
average  specific  heat  ratio  /  .  Then  the  basic  equation  for  the  shock  tube 

becomes 


A£c  is  the  heat  released  per  unit  mass  by  the  chemical  reaction  at  constant 
temperature  and  pressure  (see  Appendix  V). 

The  resulting  equations  will  be  simplified  by  using  and  Pl8  to 
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represent  shock  strengths  rather  than  M,  and  M 

si  ’# 

relations  are  applied,  it  is  found  that 


P7.'' 


"  (rJ  p„ 


When  the  normal  shock 


(38) 


and 


t 


+ 


1L. 


P  ♦  / 

J-l 


(39) 


Now  apply  a  strong  shock  approximation 
This  is  essentially  the  same  as  assuming  M  * 


Pl(  »  I  .  and  P7f  »  / 

»  I  and  M  *  |  .  Then 
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or 


=  Ff .  (*l  '  r\>  *+*  K>jr  >  P*i  *  f  >  l\i/ »  H* g) 

4ft 


(42) 


By  the  use  of  the  relation 


[  r  i  /V  l\'/x  V  5  2r  1 1/1 /nA  \H  r*~' 

&  Q  Ji  y/  (\-i )  z  p  _  <*>  /JLuJ  > 

P^'P»»Y  ^\W  (5"'  *  ^  j  vp4,/  j 


(43) 


the  dependence  of  f  upon  -~a-  could  be  replaced  by  a  dependence  upon  ^ 
only.  It  is  interesting  to  note  that  Eq.  (41)  can’t  be  appreciably  simplified 

by  making  more  stringent  assumptions  about  the  initial  shock  strength.  Of 

2 tf 

course,  it  is  true  that  for  p  ^ — - —  the  equations  are  no  longer  depend- 

ent  on  p  ,  but  this  is  the  trivial  case  corresponding  to  no  appreciable 
chemical  energy.  Suppose  that  the  initial  specific  heat  ratios  ^  ,  7^.  »  7^ 
and  the  over-all  density  ratio  are  fixed.  It  is  found  that  f^/ P^  is  still 
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a  function  of  the  five  variables  ,  <f  ,  p  ,  and  . 

To  make  this  problem  tractable  assume  that  ,  Y  ,  and  are 
known  beforehand-  Reasonable  estimates  of  these  quantities  can  be  obtained 
if  the  buffer  gas  and  its  initial  temperature  are  known.  If  the  temperature 
and  molecular  weight  of  the  driver  gas  are  known,  then 


=  K 


a  constant 


and  P  / p  depends  only  upon  p  .  The  maximizing  procedure  is  then 
straightforward  and  a  plot  of  vs.  can  be  obtained.  Of  course,  a 

different  plot  will  be  obtained  for  every  value  of  the  set  of  variables  1C  ,  7X  , 
^  ,  and  .  This  procedure  is  the  same  as  that  used  to  predict  the  first 

approximation  curves  for  an  inert-gas-buffered  tube  (Appendix  III).  By 
analogy  with  the  inert-gas  tube,  it  is  reasonable  to  assume  that  the  variation 
of  the  performance  curves  with  ^  will  be  small. 

Maximizing  (41)  for  fixed  values  of  ,  ip  ,  and  K  is  straight¬ 

forward  but  tedious.  The  results  of  the  computation  for  two  different  values 
of  K  are  shown  in  Fig.  20.  Both  curves  show  the  performance  of  a  tube  with 
a  hydrogen  driver,  stoichiometric  hydrogen- oxygen  buffer,  and  air  as  the 
driven  gas.  The  values  of  the  various  constants  used  are  as  follows: 


7i  =1.15 

=  57  kcal.  /mole 


=1.23 
y  =23.0 


for  T+  =  1000°K  X  =.22 

for  T4  -  300°K  X  =  -41 

The  X*S  represent  typical  points  from  the  "exact"  calculation  with  K  =  .22. 


The  agreement  with  the  approximate  solution  curve  indicates  that  the  assump- 
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tions  used  in  the  approximation  are  reasonable.  The  lower  curve  (  K  =  .  22) 
corresponds  to  a  heated  hydrogen  driver  and  large  values  of  for  a  given 
P  .  The  values  of  P  corresponding  to  this  curve  are  considerably 
above  the  Chapman- Jouguet  value  so  that  the  detonation  is  strongly  overdriven. 
This  performance  curve  is  almost  identical  to  that  of  a  tube  with  an  inert 
buffer  gas. 

The  curve  which  corresponds  to  a  cold  hydrogen  driver  (  K"  =  .41)  is 

somewhat  higher  than  might  be  expected.  However,  in  this  case  the  optimum 

3 

P  is  the  Chapman- Jouguet  value  up  to  about  P^  =3*10  .  This  should  give 

the  maximum  gain  over  the  inert- gas  tube  since  the  ratio  of  the  chemical  energy 

3 

to  the  total  energy  is  the  greatest  in  this  case.  Above  P4<  =3*10  the  detona¬ 

tion  is  overdriven  and  the  relative  gain  of  this  device  should  decrease.  At  very 
large  values  of  P^  the  chemical  energy  released  by  the  detonation  is  insignifi¬ 
cant  compared  to  the  total  energy,  and  the  performance  of  the  tube  is 
essentially  that  of  an  inert-gas-buffered  tube. 


3.  TESTING  TIME 


Testing  time  is  defined  as  the  interval  of  time  between  the  arrival  of  the 

shock  and  the  arrival  of  the  contact  surface  at  some  point  in  the  test  section  of 

the  shock  tube.  3y  considering  the  equation  of  continuity  we  can  arrive  at  a 
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simple  expression  for  the  ideal  testing  time. 


where  X  is  the  distance  from  the  diaphragm.  Ideally  the  testing  time  could 
be  increased  to  any  value  merely  by  building  the  driven  section  long  enough, 
course,  for  a  simple  shock  tube  with  a  driver  section  of  finite  length,  the 


(45) 


Of 
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expansion  wave  will  reflect  off  the  end  of  the  driver  and  overtake  the 
contact  surface  (see  Fig.  21).  So,  for  a  fixed  driver  length,  there  is  an 

29 

optimum  length  of  the  driven  section  which  will  yield  the  maximum  test  time. 

The  wave  diagram  analysis  for  the  buffered  shock  tube  becomes  much 
more  complex  {Fig.  22).  Many  different  wave  interactions  must  be  considered: 
for  example,  the  interaction  of  and  C,  »  the  interaction  of  C,  and  ,  the 
interaction  of  R  and  ,  etc.  Figure  22  represents  one  possible  length  dis¬ 
tribution  that  will  lead  to  relatively  long  testing  times  in  the  driven  section.  If 
area  changes  are  present  in  the  shock  tube,  the  wave  diagram  is  further  com¬ 
plicated  by  more  reflected  waves. 

Unfortunately,  the  actual  testing  time  is  not  accurately  determined  by  the 
idealized  wave  diagram  technique.  In  the  first  place,  if  the  driven  section  is 
very  long,  shockwave  attenuation  will  be  important.  The  testing  times  calculated 
for  ideal  gases  will  be  lowered  somewhat  by  considering  real  gas  effects.  How¬ 
ever,  Eq.  (45)  does  not  require  that  the  driven  gas  be  ideal,  so  that  real  gas 
calculations  can  be  used  with  it.  In  practice  the  contact  surface  will  not  be  a 
discontinuity  but  will  be  turbulent  and  diffuse.  This  reduces  the  test  time  con¬ 
siderably;  a  good  rule  of  thumb  is  to  multiply  the  ideal  test -time  by  a  factor  of 
one -half. 

At  pressures  below  about  five  millimeters  of  mercury  the  testing  time  is 

found  to  be  much  lower  than  that  predicted  even  with  the  above  corrections. 

Roshko  showed  that  this  could  be  largely  explained  in  terms  of  the  loss  of  mass 
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across  the  contact  surface  due  to  the  laminar  boundary  layer  behind  the  shock. 

He  solves  the  problem  in  terms  of  two  similarity  parameters  X  and  T 
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(46) 


X  - 


Jj_  (h,-/) 

r,  f\, 


r  =  x  (r„-i) 


r  **S, 


't 


where  S.t.p. 

O' 

r 

*2 

cL 

A 

Roshko  found 


implies  room  temperature  and  atmospheric  pressure 
parameter  depending  on  P  ,  y  ,  and  (T 
Prandtl  number 
actual  testing  time 
compressibility  of  shocked  gas 
diameter  of  driven  section 
dynamic  viscosity  coefficient 
the  relation  between  the  two  parameters  to  be 


] 


(47) 


(48) 


for  which  T  — >*  1  as  X  — >  op  .  Hooker  expanded  Roshko' s  theory  to  account 

for  boundary  layer  mass  accumulation  between  the  shock  wave  and  the  contact 
22 

surface.  Hooker's  equation  for  X  (T)  is  more  complicated  than  Eq.  (47),  but 

the  improvement  is  relatively  minor. 

Recently  Mirels  has  analyzed  this  problem  and  removed  some  inconsis- 
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tencies  which  appeared  in  the  previous  analyses.  In  particular,  by  consider¬ 
ing  the  variations  in  free  stream  conditions  between  the  shock  and  the  contact 
surface  he  was  able  to  find  improved  values  of  the  parameter  £  ,  The  im¬ 
provement  provides  much  better  correlation  with  experimental  data  for  low 
shock  Mach  numbers;  for  large  Mach  numbers  the  analysis  basically  agrees 
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with  that  of  Roshko  and  Hooker  and  with  the  experimental  data. 

The  fact  that  “f  approaches  a  limit  indicates  that  the  actual  testing 
time  does  not  continue  to  increase  as  the  length  of  the  driven  section  is  in¬ 
creased.  After  a  certain  time  the  contact  surface  moves  at  the  same  speed 
as  the  shock,  and  the  testing  time  is  at  a  constant  value  that  shall  be  denoted 
by  Tnf|>,  .  Roshko  shows  that 


_T, 


(49) 


Throughout  the  analyses  of  Roshko,  Hooker,  and  Mirels  it  is  assumed 

that  the  boundary  layer  is  thin  (relative  to  the  tube  hydraulic  radius)  and 

laminar  between  the  shock  and  the  contact  surface.  Mirels  has  evaluated  both 

the  Reynolds  number  and  the  boundary  layer  thickness  at  the  contact  surface 

(when  it  is  at  the  maximum  distance  from  the  shock)  and  found  these  assumptions 

to  be  good  for  strong  shocks  and  low  initial  pressures.  For  relatively  weak 

shocks  (under  Mach  three,  say)  and  low  initial  pressure  such  as  those  used  by 
20 

Duff,  the  boundary  layer  becomes  rather  thick  at  the  contact  surface. 

For  a  fixed  value  of  shock  Mach  number  the  maximum  testing  time  as 
determined  by  boundary  layer  development  depends  directly  upon  the  initial 
pressure  and  on  the  diameter  squared.  This  fact  has  aroused  interest  in  shock 
tubes  with  large  diameter  driven  sections  for  improving  testing  time.  The 
size  of  the  shock  tube  driver  is  usually  limited  by  structural  and  economic  con¬ 
siderations.  Therefore,  shock  tubes  with  larger  diameter  driven  sections  than 
driver  sections  have  been  studied. 

A  simple  shock  tube  with  an  increase  in  area  at  the  diaphragm  admits  to 
one  of  four  possible  flow  patterns  (see  Appendix  II  and  Fig.  23).  Generally, 
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for  strong  shocks  the  flow  will  be  of  type  (a),  but  for  very  large  area  changes 
the  other  types  are  possible.  Figure  24  illustrates  the  performance  of  such  a 
tube  with  varying  values  of  driven  pressure  and  diameter.  With  fixed  driver 
conditions  (hydrogen,  =  750°F,  -p /  =  15,000  psia,  ^  =  3.5")  and 

the  shock  Mach  number  held  at  20,  the  pressure  in  the  driven  section  must 
be  lowered  for  operation  at  larger  diameters.  These  calculations  were  performed 
assuming  that  hydrogen  is  an  ideal  gas  with  /  =7/ 5  and  using  the  properties 
of  equilibrium  air.  From  Fig.  24  it  can  be  seen  that  testing  time  is  being  lost 
by  going  to  larger  diameter  driven  sections  since  the  quantity  n**  decreases. 
The  maximum  testing  time  is  obtained  with  a  constant-area  buffered  tube. 

The  analysis  of  buffered  tubes  with  larger  driven  sections  than  driver 
sections  is  more  complex.  For  example,  if  all  the  area  change  is  at  the  second 
diaphragm,  a  typical  wave  diagram  is  shown  in  Fig.  25.  The  upstream-facing 
wave  W  is  either  a  shock  wave  or  a  rarefaction  fan  depending  upon  the  final 
shock  strength  and  the  area  ratio.  The  flow  from  region  2  to  region  5  is  a 
steady  expansion  and  is  completely  determined  by  conditions  in  region  2  and  the 
area  change.  The  most  efficient  operation  of  this  type  of  shock  tube  occurs 
when  the  buffer  pressure  is  so  chosen  that  W  reduces  to  a  sound  wave.  Of 
course,  if  the  area  change  is  small,  it  may  not  be  possible  to  reduce  it  to  a 
sound  wave.  Figure  24  compares  the  performance  of  the  shock  tube  configura¬ 
tion  of  Fig.  25  at  the  optimum  buffer  pressure  with  that  of  a  simple  shock  tube 
with  an  area  change  and  the  same  driver  conditions.  The  comparison  certainly 
indicates  the  importance  of  using  a  buffer  with  this  type  of  a  shock  tube. 

Using  the  values  of  -p>  ,  >  and  M.  from  Fig.  24,calculations  of 

were  made  from  Roshko's  theory.  The  testing  time  drops  with 
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increasing  diameter,  though  not  as  severely  as  for  the  simple  shock  tube. 

For  certain  arc-heated  facilities  the  testing  time  is  found  to  actually  increase 

5 

with  increasing  diameter  and  fixed  Mach  number.  With  an  arc-heatcd  helium 
driver  the  speed-of- sound  ratio  from  the  driver  gas  to  the  driven  gas  may  be 
as  high  as  thirteen.  In  this  range  the  gain  in  shock  Mach  number  with  in¬ 
creasing  pressure  ratio  is  greater  than  in  the  more  conventional  range  of  speed- 
o£-sound  ratios.  Nevertheless,  the  gain  in  testing  time  is  quite  small.  There 
appears  to  be  no  strong  reason  for  using  large  diameter  driven  sections  when 
operating  at  high  Mach  numbers. 

Now  consider  a  buffered  tube  similar  to  that  of  Fig.  25  but  with  a 

converging-diverging  area  change  at  the  second  diaphragm.  A  typical  wave 

diagram  for  this  configuration  is  shown  in  Fig.  26.  Again  W  may  be  a  shock 

or  an  expansion,  and  this  must  be  determined  in  each  particular  problem. 

Figure  27  compares  the  performance  of  such  a  configuration  with  that  of  a  tube 

with  just  a  monotonic  divergence  at  the  second  diaphragm.  Notice  that  there 

is  always  a  loss  in  performance  by  using  the  converging-diverging  section. 

This  result  is  analogous  to  the  result  obtained  for  the  case  of  an  area  change 
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at  the  first  diaphragm  (see  Appendix  II).  It  is  apparently  true  that  the  most 
efficient  way  to  accomplish  an  area  change  at  a  diaphragm  is  monotonically. 

Finally,  examine  the  effect  of  distributing  a  total  area  expansion  between 
the  two  diaphragm  stations.  Figure  28  presents  results  for  given  over-all 
conditions  and  varying  buffer  diameter  and  pressure.  The  maximum  final 
shock  Mach  number  increases  as  the  buffer  diameter  is -increased  from  the 
driver  size  to  the  driven  section  size.  Although  the  performance  continues 
to  improve  as  the  buffer  diameter  is  further  increased,  the  improvement  is 
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rather  limited  and  by  the  time  that  A,  =  2 A.  the  performance  curve  is  again 
below  the  curve  for  At~  Ag  •  So,  for  practical  purposes,  the  most  efficient 
geometrical  configuration  is  one  for  which  all  the  area  expansion  occurs  at 
the  first  diaphragm. 


4.  CONCLUDING  REMARKS 

A  method  has  been  presented  for  quickly  determining  the  optimum  operating 
conditions  of  buffered  shock  tubes  and  the  final  shock  strength  obtainable  at  these 
optimum  conditions.  While  the  resulting  equations  are  based  on  a  strong  shock 
assumption,  these  equations  are  quite  accurate  for  most  of  the  cases  of  interest. 
Even  in  problems  to  which  the  assumptions  do  not  wholly  apply,  the  analysis 
indicates  the  range  of  conditions  to  consider  for  making  exact  calculations. 

Most  of  the  report  is  concerned  with  applying  the  strong  shock  analysis 
to  different  buffered  shock  tube  configurations.  Although  much  has  been  written 
about  buffered  shock  tubes,  it  seems  worthwhile  to  summarize  our  present 
knowledge  about  these  devices. 

1.  For  given  over-all  conditions  a  buffered  shock  tube  will  produce 

stronger  final  shocks  than  a  simple  shock  tube  provided  that  the  optimum 
buffer  conditions  are  used.  However,  if  the  over-all  density  ratio  ts  of  order 
ten,  the  improvement  is  negligible.  The  relative  improvement  increases  as 
the  over-all  density  ratio  approaches  its  two  limits  (r4,->o,  ■ 

2.  Within  the  framework  of  the  strong  shock  approximation,  the  final 
shock  strength  2_  depends  only  upon  the  upstream  density  ratio  p  when  the 
specific  heat  ratios  and  the  over -all  density  ratio  are  specified.  After  optimiz¬ 
ing,  2  -maximum  depends  only  upon  P  and  the  specific  heat  ratios. 

3.  Other  factors  being  equal,  the  best  driver  gas  is  diatomic  and  the 
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best  buffer  gas  is  monatomic.  The  specific  heat  ratio  of  the  driven  gas  is 
relatively  unimportant. 

4.  If  the  strong  shock  assumption  is  relaxed  somewhat,  it  is  found  that 

2  -maximum  depends  not  only  upon  but  also  (to  a  lesser  degree)  upon 

fp  Y/l  ’ 

K  =  I-—*1  ■}  .  For  fixed  P  ,  j*  -maximum  increases  as  K  increases 

\P4,  /  4ft  a 

(  decreases). 

5.  The  performance  of  a  buffered  shock  tube  improves  as  the  number  of 
buffers  increases  provided,  of  course,  that  the  optimum  conditions  are  used. 
However,  as  the  number  of  buffers  increases,  the  final  shock  strength  approaches 
its  asymptotic  value,  and  the  relative  gain  decreases. 

6.  A  buffered  tube  with  an  area  contraction  at  either  or  both  diaphragms 
yields  higher  shock  strengths  than  a  constant-area  tube.  The  optimum  distri¬ 
bution  of  a  given  area  contraction  is  in  equal  proportions  between  the  two 
diaphragm  stations.  However,  if  all  the  area  reduction  is  to  be  at  one  diaphragm 
station,  it  is  better  for  the  change  to  be  at  the  second  diaphragm. 

7.  A  shock  tube  with  a  detonable  gas  in  the  buffer  can  produce  stronger 
final  shocks  than  inert-gas-buffered  tubes.  The  maximum  gain  is  obtained 
when  the  initial  shock  (detonation)  is  near  the  Chapman- Jouguet  value.  However, 
for  the  production  of  very  strong  final  shocks,  the  detonation  must  be  consider¬ 
ably  overdriven.  In  these  cases  the  gain  in.  performance  over  the  inert-gas 
tubes  is  not  significant. 

8.  For  fixed  driver  conditions,  increasing  the  diameter  of  the  driven 
section  requires  that  the  pressure  in  the  driven  section  be  lowered  in  order  to 
maintain  the  same  final  shock  Mach  number.  The  net  effect  of  the  increased 
diameter  but  decreased  pressure  is  to  reduce  the  available  testing  time  in 
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most  cases  of  interest.  When  the  speed-of- sound  ratio  across  the  diaphragm 
is  very  high,  only  a  small  increase  in  testing  time  may  be  obtained. 
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APPENDIX  I 


DERIVATION  OF  THE  BASIC  EQUATIONS  FOR  CONSTANT - 
AREA  BUFFERED  SHOCK  TUBES 

Many  of  the  results  derived  in  these  appendices  are  well  known  and  are 
presented  in  several  places  in  the  literature.  They  are  included  here  for  con¬ 
venience  and  to  introduce  the  notation. 

Consider  the  constant-area  buffered  shock  tube  shown  in  Fig.  1.  The  wave 
diagram  is  drawn  subject  to  the  following  assumptions:  the  flow  is  one¬ 
dimensional  and  inviscid  with  no  heat  conduction,  the  diaphragms  _re  instantan¬ 
eously  removed,  no  mixing  occurs  at  the  contact  surfaces,  and  there  is  no  shock 
wave  attenuation.  Across  the  contact  surfaces  velocities  and  pressures  are 
equal,  while  other  thermodynamic  variables  may  be  discontinuous. 

'P'x,  -  1*$  i  ft*?- 7  and  U6  =  U7  (A-l) 

The  ordinary  Rankine-Hugoniot  shock  relations  hold  across,  the  shock  waves  in 
coordinate  systems  at  rest  with  respect  to  the  shocks.  For  the  initial  shock 

p,  us,  =  (u,t  -  «, )  <A-2> 


Pi  +  P>  US,  =  7a’*.+A  K  - 


(A  -  3 ) 


+4  u.1  =-*-.  +  ^r-  +j  &,~US 


'  P,  2 


P* 


(A-4) 


where  JL  is  the  internal  energy  per  unit  mass.  A  similar  set  of  equations  holds 


for  the  final  shock  S#  .  Applying  the  method  of  characteristics  to  the  one- 

^ — 

dimensional  unsteady  motion  in  the  upstream-facing  expansion  waves  R(  and 
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j^T  it  is  found  that  the  P  Riemann  invariant  U  +  is  constant.  ^ 

Given  the  initial  conditions  in  regions  4,  l.and  8  and  the  thermal  and  caloric 
equations  of  state  for  the  gases,  the  shock  Mach  numbers  and  M,  can  be 
uniquely  determined.  In  what  follows  it  will  be  assumed  that  the  gases  used 
in  the  shock  tube  are  ideal  and  have  constant  specific  heats.  In  this  case,  the 


Riemann  invariant  expression  can  be  integrated  to  yield 

2a. 

U  +  -  =  const. 

r-t 


(A-5) 


across  an  upstream  facing  wave  (Q-wave).  Also,  since  the  flow  across  any 
expansion  is  isentropic, 


5-  •  fef 


(A-6) 


I’he  velocity  U\  is  zero  so  that  Eq  (A-5)  applied  to  R.  yields 


&+  -  + 


v-/ 


(A-7) 


From  Eqs.  (A-l),  (A-6),  and  (A-7) 


=  *  L  ‘  ~  -sg 


The  pressure  ratio  and  the  density  ratio  across  the  initial  shock  are 


-(r,-D 


(A-8) 


(A-9) 


fc+0Mti 

rz,  "  (rri)nsyz 


(A  - 1  0) 
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Using  Eq.  (A-10)  in  (A-2),  it  is  seen  that 


“2  = 


2*,,  ms;  -/ 


*i*l  -M,, 


Finally,  using  (A-8),  (A-9),  and  (A-ll),  the  well-known  simple  shock  tube 
equation  is  obtained. 


p*,= 


r,  +  i 


I- 


n  \-i 


*,*•  Ws, 


Across  the  following  fundamental  relations  hold: 


(A-12) 


7,-/  /,-/ 

=  **  +~T~  “<> 


(A-13) 


and 


2f, 

'*4  \ 

P~  =  1— 


(A-14) 


Relations  similar  to  (A-9)  and  (A-ll)  hold  across  the  final  shock  S 


s 


2  ^  M*.  ~(ra~^ 


i's  +  I 


(A-l  5) 


2  a 


U7  * 
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M  *  -  / 
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(A- 16) 


Using  (A-13),  (A-ll)  and  (A-l6)  in  (A-14),  one  finds  that 


P  = 

2-7 
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(A-17) 
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From  (A-15),  (A-17),  (A-12)  and  (A-9)  the  fundamental  equation  for  the 
constant-area  buffered  shock  txibe  is  obtained. 


p  =  P  p  p  - 
48  r7a  \i  rn 


2r,  M  *  -  (y,-;) 


y»+' 


Vj  Mi,-; 

V'  A„MSl 


r  ««  a,  ;  \ 

l  ai  *>.,  / 


'  77  K'ijf 


With  the  use  of  the  normal  shock  relation 


*•  "  ft")  M. 


and  the  ideal  gas  relation 


*  V'7^  f 


(A -18)  can  be  written  as  follows: 
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Next  consider  a  double-buffered  tube  (of  the  "unsteady- expansion'1  type) 


as  shown  in  Fig.  11.  The  wave  diagram  represents  only  the  initial  flow 
phenomena  in  such  a  tube.  Interactions  between  C  and  or  between  Cz 
and  will  generate  waves  that  will  limit  the  testing  time  (i.  e. ,  the  time 


between  the  passage  of  &  and  the  arrival  of  the  contact  surface 


The  flow  between  regions  7  and  10  is  completely  similar  to  the  flow 
between  regions  2  and  6,  so  that,  in  analogy  to  Eq.  (A-17), 
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(A-22) 


Also,  the  basic  shock  relations  across  Sja  are  similar  to  Eqs.  (A-9),  (A-ll), 
and  (A-i'9).  Therefore,  the  basic  equation  for  the  double-buffered  shock  tube 
can  be  written  down. 
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Of  course,  the  simple  shock  tube  Eq,  (A- 12)  and  the  constant-area 
buffer  Eq.  (A-21)  are  also  valid  for  this  configuration,  so  that  there  are 
three  equations  to  determine  M  ,  M  ,  and  M,  •  The  above  procedure  can 

3i  *8 

be  extended  to  a  constant-area  shock  tube  with  a  large'  number  of  buffer  sec¬ 
tions.  If  the  final  driven  section  is  labeled  fit  A-  .  the  general  equation 
(analogous  to  (A -23)  )  can  be  written  down  immediately. 


APPENDIX  II 


THE  SIMPLE  SHOCK  TUBE  WITH  AREA  CHANGE  NEAR  THE  DIAPHRAGM 

While  the  flow  in  a  duct  of  varying  cross-sectional  area  is  not  strictly 

one -dimensional,  the  simplification  of  quasi-one-dimensional  flow  in  which 

the  area  A  (x,  t)  appears  as  one  of  the  variables  in  the  fundamental  equations 

is  often  made.  Once  A  (x,  t)  is  specified,  the  problem  can  be  solved  by  the 

30 

application  of  well-developed  wave  diagram  techniques.  A  further  approx¬ 
imation  which  can  often  be  made  is  to  replace  the  gradual  change  in  duct  area 
by  a  discontinuous  change.  Such  an  approximation  can  only  be  used  to 
evaluate  the  effect  of  the  area  change  at  a  sufficient  distance  from  the 
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31  -33 

change.  Since  the  changes  in  the  flow  variables  are  assumed  to  occur 

instantaneously,  they  may  be  computed  by  using  the  steady  flow  equations  for 
an  isentropic  nozzle. 

The  possible  wave  patterns  for  a  simple  shock  tube  with  an  area  change 
near  the  diaphragm  that  are  compatible  with  the  above  assumptions  are  shown 
in  Fig.  23.  The  wave  diagrams  always  indicate  a  primary  shock,  an  upstream¬ 
facing  expansion  wave,  and  a  contact  eurface.  In  general,  an  upstream-facing 
secondary  wave  will  be  required  to  match  conditions  at  the  interface.  The 
steady  flow  relations  are  assumed  to  hold  from  region  3a  to_3b.  In  any  partic¬ 
ular  case  the  flow  pattern  will  be  determined  by  the  initial  conditions  and  the 
area  ratios. 

Russell  has  done  a  detailed  study  of  all  four  types  of  solutions  for  varying 
values  of  A A,  and  A*/^,  .  33  He  demonstrates  that  for  any  value  of  A^/a%  , 

the  shock  strength  increases  with  increasing  A4/A,  ,  and  for  any  value  of 
A4/A,  the  maximum  shock  strength  corresponds  to  the  maximum • 
Therefore,  the  best  over -all  performance  is  obtained  with  a  monotonic  area 
contraction  at  the  diaphragm  station.  With  this  configuration  only  solutions  of 
type  (a)  or  type  (d)  are  possible.  For  these  solutions  the  steady  expansion  is 
subsonic  and  the  second  non- steady  expansion  (if  any)  is  supersonic.  The  im¬ 
proved  performance  follows  from  the  fact  that  a  steady  subsonic  expansion  is 

more  efficient  than  an  unsteady  subsonic  expansion  for  the  conversion  of 

31  32 

thermal  to  kinetic  energy.  ’ 

The  following  analysis  is  valid  for  any  flow  of  either  type  (a)  or  type  (d) 
whether  Af-  A#  or  not.  For  the  steady  flow  from  region  3a  to  3b  the 
following  relations  hold: 


49 


(A-25) 


(A-26) 


(A-27) 


(A-28) 


(A -29) 
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Define 
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Equation  (A -31)  can  be  written  as 
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Equations  (A -9)  and  (A-ll)  may  be  used  to  obtain 
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(A-30) 


(A-31 ) 
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Then  Eq.  (A- 34)  is  merely  the  simple  shock  tube  Eq.  (A -12)  for  a  shock  tube  with 
a  pressure  ratio  and  a  speed  of  sound  ratio  . 

The  specific  case  of  interest  is  the  configuration  for  which  ^  . 

Then  if  ^  j  ,  M,b=  /  “and  <£  can  be  calculated  from  (A-26)  and  (A-30), 

If,  on  the  other  hand,  M.  I  ,  M,.  =  M,  and  the  problem  must  be  solved  by 
iteration. 

19 

Hall  and  Russo  have  pointed  out  that  for  Ms  exceeding  about  three,  and 
for  a  given  driver  gas  specific  heat  ratio,  shock  tubj  performance  can  be  repre¬ 
sented  by  a  single  curve  for  all  initial  conditions.  If  Ms  is  greater  than  about 
three,  certain  terms  in  Eq.  (A-34)  can  be  neglected.  Then  the  equation  can  be 


Define  two  new  variables  by 


(A-36) 


(A-37) 


The  basic  equation  has  been  reduced  from  a  relation  involving  six  variables  M,(  , 
P41  ,  ^  ,  P  ,  7,  ,  and  to  a  relation  among  four  variables  "if ,  P^  , 

7,  ,  and  / ^  .  Figure  6  shows  a  plot  of  vs.  for  various  values  of  driver 


gas  specific  heat  ratio.  The  use  of  different  values  for  7,  does  not  noticeably 


52 


affect  the  curves  of  Fig.  6,  so,  the  variation  in  performance  with  driven  gas 
specific  heat  ratio  can  be  neglected.  Figure  6  indicates  that  for  given  K  >  2, 
increases  as  decreases.  By  using  numerical  values  it  is  found  that  for 
given  P  ,  (the  actual  shock  strength  parameter)  increases  as  7^ 

decreases. 


APPENDIX  III 


CALCULATIONS  OF  THE  OPTIMUM  PERFORMANCE  OF 
CONSTANT -AREA  BUFFERED  TUBES 

The  performance  of  a  constant-area  buffered  tube  can  be  calculated  from 
the  buffer  Eq.  (A-21)  and  the  simple  shock  tube  Eq.  (A -12).  In  the  main  body 

X  X 

of  the  text  it  is  shown  that  when  M  and  M  are  large  compared  to  one, 

si  5s 

approximation  1  is  obtained. 
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The  additional  assumption  that  , 

5 ,  “  / 


(A-38) 


simplifies  (A-38)  to  approx¬ 


imation  2,  which  no  longer  explicitly  depends  upon  M3 

O  y  1 

be  n'A i  W 

,*•  -  -  *  )  _  o  c.,P4ktt 


2:  =  C.H 
a  i  / 


-ea  <r 


(A-39) 


2  can  be  maximized  with  respect  to  2  by  differentiating  (A-39)  with  respect 

8  / 


to  2  and  setting 


2L 

5>iL 


equal  to  zero.  Upon  differentiation,  Eq.  (A-39) 
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becomes 
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The  possibility  that 
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BE,  yt 
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(A -41) 


J,+l  ' 

is  excluded  because  this  would  imply,  by  (A-39),  that  =  0  or  p  becomes 

infinite.  Therefore,  (A-40)  can  be  simplified  to 

BE^g, _  _ 
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(A- 43) 
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where 


The  use  of  fA-431  to  eliminate  7  from  (A -391  leads  to 


This  can  be  solved  for  . 


(A -44) 


(A -45) 


(A-46) 


With  (A-46)  and  (A-43)  one  finds  that 


_  >A 

*8  ~  +X, / 


(A-47) 


The  possibility  of  determining  the  maximum  value  of  j*  from  approxi- 
mation  1  now  will  be  considered.  Equation  (A -38)  can  b^r  written  symbolically 


z»  =  f5  fr' V Vz.’MvrJ 


(A-48) 


So,  lor  fixed  specific  heat  ratios  and  over-all  density  ratio  the  conditions  for 
the  maximum  are 
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—  0  and 


=  o 


(A -49) 


Equation  (A-38)  can  be  written  as 


=  c,  ^<1*  f  (M  J  - 


where 


f/M  N 

Differentiating  (A-50)  with  respect  to  2,  leads  to 
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and  differentiating  with  respect  to  Ms  leads  to 
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Equations  (A-50),  (A-51)  and  (A-52)  can  be  combined  to  yield 


(A-50) 


(A-51) 


(A-52) 
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V/hen  values  of  are  used  to  evaluate  these  expressions,  it  is  found 

that  the  only  possible  values  of  are  those  in  the  range  of  about  2  to  17. 
These  are  the  only  values  of  P^$  for  which  the  expression  in  brackets  in 

Eq.  (A-50)  can  be  equal  to  one.  For  larger  values  of  p  (which  are  usually 

48 

the  ones  of  interest)  the  maximum  value  of  £  cannot  be  determined  by  the 

8 

above  procedure. 

A  procedure  to  determine  the  optimum  conditions  for  approximation  1 

curves  such  as  appear  in  Figs.  10  and  11  will  now  be  considered.  These  curves 

were  drawn  subject  to  the  constraint  that  the  temperatures  and  molecular 

weights  are  specified,  and,  therefore,  £  is  no  longer  a  function  of  two 

8 

variables.  The  equation  for  approximation  1  may  be  expressed  as  follows: 
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The  constraint  can  be  expressed  as 


£41  _  yC*-+!  _  ^3. 
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Therefore, 
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The  optimum  condition  can  be  found  from  the  equation 

3Ft  _ 


(A-59) 


(A-60) 


The  calculation  is  straightforward  but  quite  involved,  so  it  will  not  be  presented. 
The  most  important  result  of  the  above  analysis  is  the  following:  for 

given  p  >  ?.  -maximum  increases  as  J<  increases  (  M  decreases). 

4#  c  *1 


APPENDIX  IV 

THE  BUFFERED  SHOCK  TUBE  WITH  AN  AREA  CONTRACTION 
AT  THE  SECOND  DIAPHRAGM 


The  use  of  an  area  contraction  in  conjunction  with  a  buffered  shock  tube 
has  been  considered  by  several  authors.  *  ’  ®  If  the  area  change 

is  located  at  the  first  diaphragm,  the  flow  is  that  of  a  simple  shock  tube  with 
an  area  change  followed  by  a  flow  identical  to  that  in  a  constant-area  buffered 
tube.  This  case  is  fully  considered  in  the  main  body  of  the  report  (Sec.  2.2.  3). 

Suppose  that  an  area  contraction  exists  at  the  second  diaphragm  station 
of  a  buffered  shock  tube  and  that  the  diaphragm  ruptures  immediately  upon 
impact  of  the  initial  shock  (Fig.  16).  If  the  contraction  is  approximated  by  a 
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discontinuous  change,  the  resulting  simplified  wave  diagram  will  include  the 
incident  and  transmitted  shocks,  an  interface  that  separates  the  two  gases, 
and  a  reflected  wa/e.  The  cases  of  interest  will  involve  supersonic  flow  behind 
the  initial  shock  (m  *>i)  ,  and  under  this  assumption  three  possible  solutions 
exist.  If  the  contraction  is  sufficiently  small,  the  reflected  wave  is  an  upstream  - 
facing  expansion  (Q-wave)  which  is  swept  downstream  (Fig.  16-a).  At  larger 
contraction  ratios  the  flow  is  reduced  to  the  sonic  value  and  a  shockwave  is 


required  to  match  conditions  at  the  interface.  For  a  limited  range  of  contrac¬ 
tion  ratios  the  shock  can  be  stationary  at  the  area  change  (Fig.  16-b),  but  as 
A,  j  Ag  increases  a  reflected  shock  is  found  (Fig.  16-c). 

There  has  been  a  considerable  amount  of  interest  in  this  problem  in 
34  35 

recent  years.  ’  For  certain  values  of  incident  shock  strength  and  contrac¬ 
tion  ratio, all  three  of  the  above  solutions  will  satisfy  the  conditions  of  the  pro¬ 
blem.  This  "region  of  ambiguity"  was  found  to  be  a  direct  result  of  the  simplify¬ 
ing  assumptions  implicit  in  the  wave  diagram  technique.  More  exact  analyses 
indicate  which  solution  will  hold  in  any  particular  case. 

For  cases  of  practical  interest  the  area  change  will  be  so  large  that  a 
reflected  shock  exists,  and  the  flow  is  adequately  represented  by  Fig.  16-c. 

The  flow  from  region  4  to  region  1  will  obviously  be  represented  by  the  simple 
shock  tube  equation.  After  the  passage  of  the  reflected  shock  the  flow  will 
have  a  small  subsonic  velocity  U  5  •  In  a  coordinate  system  which  is  station¬ 
ary  with  respect  to  the  reflected  shock 


(A-61) 
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where 


speed  of  reflected  shock 


reflected  shock  Mach  number 


The  ordinary  shock  relations  hold  across  the  reflected  shock. 

2*/  -  fc-f ) 


(A- 62) 


0?  r,  -  fc  -  /)]'A  [fc  -  0  Msl  -6  2]  'A 

(V'>  Ms 


(A-63) 


The  flow  from  region  5  to  the  beginning  of  the  unsteady  expansion  51  is  steady. 

For  a  monotonic  convergence  =1,  but  if  the  area  change  is  a  converging- 

/ 

diverging  section  M,  >  1  .  The  usual  steady  flow  relations  apply. 
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Define  &  by  the  expression 
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Applying  the  isentropic  relation. 
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it  is  found  that 
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(A-69) 


(A-70) 


(A-71) 


Equation  (A-71)  may  be  interpreted  as  meaning  that  the  flow  from  region  5  to 

region  8  is  equivalent  to  the  flow  in  a  simple  shock  tube  with  a  pressure  ratio 

r,-l 

GP.  and  a  speed  of  sound  ratio  G  '  .  Therefore,  G  is  an 

vU  *  ^  O 

equivalence  factor  for  an  area  change  at  the  second  diaphragm.  It  is  import¬ 
ant  to  note  that  this,  is  not  the  same  as  the  "  "  factor  used  previously.  If 

equations  (A-62),  (A-63),  and  (A-71)  are  combined,  the  basic  equation  for 
this  type  of  shock  tube  is  obtained. 
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However,  the  equation  still  involves  two  unknowns,  and 

The  equation  of  continuity  across  the  reflected  shock  is 


or 


'“l  (Uxr+Ul)  =  <°5  (U5r  +  Us) 


By  using  the  relation 
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in  (A-74),  it  is  found  that 
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M  is  known  and  can  be  calculated  from  Eq.  (A-65),  and  if  (A-65)  is 
used,  (A-76)  becomes  a  fourth  order  algebraic  equation  for  Ms  •  In  any 


(A-76) 
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practical  case  Eq.  (A-76)  can  be  solved  by  iteration. 


If  the  area  contraction  is  very  large,  Us  ^ 
reflected. 
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2r,  M,,  -  fa-/) 


This  expression  may  be  used  to  eliminate 
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0  and  the  shock  becomes  fully 

(A-77) 

from  Eq.  (A-72). 

X 


where 


since  Ms  =  0  and 


Mr  =1  for  an  infinite  monotonic  convergence. 


APPENDIX  V 

DERIVATION  OF  THE  DETONATION  BUFFER  EQUATIONS 

Under  normal  conditions  a  flame  in  a  tube  filled  with  combustible  gas 
will  propagate  at  the  low  velocity  of  a  few  meters  per  second.  Under  certain 
conditions,  however,  this  slow  combustion  process  is  changed  into  a  very 
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rapid  process  that  propagates  at  supersonic  speeds.  This  second  type  of 
combustion  process  is  known  as  a  detonation.  In  the  following  it  will  be 
assumed  that  the  chemical  reaction  occurs  instantaneously  across  a  sharply 
defined  front.  Then  the  only  difference  between  a  detonation  and  a  pure  shock 
wave  is  that  the  chemical  nature  of  the  burnt  gas  differs  from  that  of  the  un¬ 
burnt  gas  and  that  the  reaction  influences  the  energy  balance. 

In  the  following  three  pages  the  fundamentals  of  detonation  wave  theory 

25  28 

will  be  reviewed.  Similar  developments  can  be  found  in  standard  texts.  ’ 
Referring  to  a  detonation  moving  into  a  quite  gas  at  speed  {J  ,  the 


Rankine-Hugoniot  relations  become 
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AE,  is  the  heat  of  chemical  reaction  per  unit  mass  at  constant  temperature 

€ 

and  pressure.  From  (A-79)  and  (A-80)  the  relation 


-Px-f'i 


P,  V,  =  Pt  v2  >0 


(A-8Z) 


follows,  just  as  in  the  case  of  pure  shock  waves.  Here  /ir  =-1/^°  is  specific 
volume.  Equation  (A-82)  implies  that  there  are  two  different  types  of  processes 
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compatible  with  the  conservation  laws:  those  in  which  both  pressure  and 
density  increase  and  those  in  which  both  pressure  and  density  decrease. 
Processes  of  the  first  kind  are  detonations  while  processes  of  the  second 
kind  are  slow  combustions  or  deflagrations.  By  considering  pure  shockwaves 
(A£c-0)  ,  it  is  found  that  deflagration-type  processes  ^t>^are  excluded 

because,  they  would  involve  a  decrease  in  entropy.  However,  this  argument  does 
not  apply  to  the  case  where  a  chemical  reaction  occurs. 

By  proper  manipulation  of  (A -79)  and  (A- 80)  the  relation 


\  -  v, 2  =  (i>r  px)  fc  +•  'vj 


(A-83) 


is  found,  which  when  used  in  Eq.  (A- 81)  yields 

■s.- =x  (px  +  P')(nr,-/n)  +  AEe  (A-84) 

This  relation  involves  only  thermodynamic  variables  and  is  commonly 
called  the  Hugoniot  relation.  The  internal  energy  can  be  considered  a  function 
of  the  pressure  and  specific  volume;  in  particular  for  an  ideal  gas  with  constant 
specific  heats 

X  "  r-\  (A-85) 


Therefore,  knowing  ^.Ec,  Jb  •-f>i  an^  4^,  >  the  Hugoniot  relation  can  be 
plotted  in  a -^-/irplane.  This  plot  of  (A-84)  subject  to  the  restriction  of  Eq. 
(A-82)  is  called  a  Hugoniot  curve  (Fig.  29).  The  curve  has  a  detonation  branch 
(/tr£rf)  and  a  deflagration  branch  -p>t)  .  The  present  study  will  just  consider 
the  detonation  branch;  most  phenomena  to  be  considered  for  this  branch  will 
have  an  analogue  in  the  deflagration  branch. 

A  straight  line  through  the  point  will  in  general  intersect  the 

Hugoniot  curve  in  two  points.  By  increasing  the  slope  of  this  line  (decreasing 
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the  angle  $  }  the  two  points  of  intersection  eventually  come  together  at 
point  C  which  is  known  as  the  Chapman- Jouguet  point.  Detonations  repre¬ 
sented  by  points  below  C  will  be  called  weak  detonations;  those  represented 
by  points  above  C  will  be  called  strong  or  over -driven  detonations. 

The  Chapman- Jouguet  process  has  certain  peculiar  properties  that  can 
be  derived  from  thermodynamic  considerations.  Among  these  are  the  following; 
of  all  detonation  processes  the  detonation  velocity  U ^  and  the  entropy  of  the 
burnt  gas  Sx  are  minima  for  a  Chapman- Jouguet  process,  and  a  Chapman- 
Jouguet  front  when  observed  from  the  burnt  gas  moves  with  the  sound  speed. 

Vx  =  Ux  -  <tz  (A-36) 

Also,  by  thermodynamic  arguments,  certain  general  statements  can  be 
derived  which  are  called  Jouguet’s  Rules.  These  can  be  stated  as:  the  gas 
flow  relative  to  the  reaction  front  is  supersonic  ahead  of  any  detonation,  super¬ 
sonic  behind  a  weak  detonation,  and  subsonic  behind  a  strong  detonation.  If  a 
more  detailed  analysis  is  made  of  the  reaction  process,  taking  heat  conduction 
and  viscosity  into  account,  it  can  be  shown  that  weak  detonations  are  possible 
only  under  extreme  and  rare  circumstances. 

Consider  a  detonation  front  which  is  being  followed  by  a  piston  moving 
with  constant  velocity  14^,  .  This  problem  is  quite  similar  to  the  problem 

of  a  detonation  in  a  shock  tube  where  the  condition  of  constant  pressure  and 
velocity  at  the  contact  surface  replaces  the  piston.  Suppose  that  —  ax  ~ 

Then  the  velocity  of  the  gas  behind  the  front  is  greater  than  that  which  exists 
at  the  Chapman- Jouguet  point,  and  the  detonation  is  over-driven.  If  14^,=  * 

a  pure  Chapman-Jouguet  detonation  is  sustained.  If  U.^4.  —  >  a 

Chapman -Jouguet  detonation  will  occur,  but  is  will  be  immediately  followed 
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by  a  centered  expansion  wave.  The  front  of  the  expansion  will  just  coincide 
with  the  detonation  front  since 


(A -87) 


If  the  piston  is  not  moved  at  all,  =  0,  and  a  Chapman- Jouguet  detona¬ 
tion  is  still  possible.  In  this  case  the  expansion  will  decelerate  the  flow  to  zero 
velocity.  That  this  detonation  is  the  one  that  actually  occurs,  is  the  Chapman- 
Jouguet  hypothesis. 

During  the  use  of  a  detonation  in  a  shock  tube  it  is  desirable  to  eliminate 
the  expansion  fan  which  follows  the  front.  Therefore,  only  driver  pressures 
which  are  great  enough  to  support  or  overdrive  a  pure  Chapman- Jouguet  deton¬ 
ation  will  be  considered.  F?om  (A -79) 
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After  multiplying  these  two  relations  together  and  dividing  by  d,  , 
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Assume  that  both  the  burnt  and  unburnt  gases  are  ideal  gases  with  constant 
specific  heats. 
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When  (A-92)  is  used  in  Eq.  (A-81), 
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From  Eqs.  (A-79),  (A-88),  (A-93),  and  the  definition  of  the  speed  of 


sound 
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If  this  expression  is  combined  with  (A-89), 
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This  can  be  used  in  (A- 94)  to  obtain  as  a  function  of  P^  ,  but  the  corres¬ 
ponding  expression  is  quite  complex.  Consequently,  we  shall  use  as  a 

measure  of  shock  strength  rather  than  M  - 

*i 

Consider  a  constant-area  buffered  tube  similar  to  that  shown  in  Fig.  1 
with  the  exception  that  the  buffer  gas  is  a  detonable  mixture.  Assume  that  the 
driver  and  driven  gases  are  perfect  gases.  The  equations  which  hold  across  the 
detonation  front  have  been  derived  above.  The  burnt  gas  will  not  behave  as  a 
perfect  gas  in  the  expansion  from  region  2  to  region  6,  and  the  flow  through 
the  fan  should  properly  be  evaluated  by  numerically  integrating  the  equations  of 
characteristics.  For  the  purpose  of  our  simplified  analysis  assume  that  the 
isentropic  relation  can  be  written  as 
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(A-96) 


r2-  =  M 

■p-,.  w  > 

where  /  is  the  average  specific  ratio  between  region  2  and  region  6,  and 

•*»r 

the  nonsteady  expansion  relations  can  be  adequately  represented  by 


(A-97j 


Then  the  basic  equation  for  this  shock  tube  is 


X 


(A-98) 
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Figure  1  INITIAL  FLOW  IN  CONSTANT-AREA  BUFFERED  SHOCK 
TUBE  -  UNSTEADY-EXPANSION  TYPE 
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PERFORMANCE  OF  CONSTANT-AREA  BUFFERED  TUBE  -  UNSTEADY- EXPANS  I  ON  TYPE 
WITH  THE  SAME  DIATOMIC  GAS  USED  THROUGHOUT  (Ref.  9) 


PERFORMANCE  OF  CONSTANT-AREA  BUFFERED  TUBE  -  UNSTEADY-EXPANSION  TYPE 
WITH  THE  SAME  MONATOMIC  GAS  USED  THROUGHOUT 
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Figure  5  EFFECT  OF  PRESSURE  RATIO  DISTRIBUTION  AND  BUFFER 
GAS  MOLECULAR  WEIGHT  UPON  STRENGTH  OF  FINAL  SHOCK 


CE  OF 
£  FOR 


DENSITY  RATIO  WITH  CONSTANT-AREA  BUFFERED  TUBE 


Figure  8  OPTIMUM  PERFORMANCE  OF  A  CONSTANT-AREA  BUFFERED  SHOCK  TUBE 
COMPARED  WITH  THE  PERFORMANCE  OF  A  SIMPLE  SHOCK  TUBE 
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Figure  9  COMPARISON  OF  THE  EXACT  EQUATION  WITH  APPROXIMATIONS 
ONE  AND  TWO  FOR  THE  CONSTANT-AREA  IUFFERED  TUBE 
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Figure  12  OPTIMUM  PERFORMANCE  OF  BUFFERED  SHOCK  TUBE  WITH  ONE  OR  TWO 
BUFFERS  COMPARED  WITH  SIMPLE  SHOCK  TUBE 
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BUFFERS  COMPARED  WITH  SIMPLE  SHOCK  TUBE 
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Figure  14  OPTIMUM  PERFORMANCE  OF  BUFFERED  SHOCK  TUBE  WITH  ONE  OR  TWO 
BUFFERS  COMPARED  WITH  SIMPLE  SHOCK  TUBE 


Figure  15  WAVE  DIAGRAM  FOR  BUFFERED  SHOCK  TUBE  WITH  AREA 
COHTRACTIOH  AT  FIRST  DIAPHRAGM  STATION 
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POSSIBLE  WAVE  PATTERN 
SECOND  D I APHRAGM  IN  A 
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PERFORMANCE  OF  VARIOUS  BUFFERED  SHOCK  TUBE  CONFIGURATIONS 
WITH  AN  OVER-ALL  AREA  CONTRACTION  OF  9.0 
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Figure  18  COMPARISON  OF  THE  EXACT  AND  APPROXIMATE  EQUATIONS  FOR  THE  BUFFERED  SHOCK  TUBE 
WITH  A  LARGE  AREA  CONTRACT I  OH  AT  THE  SECOND  DIAPHRAGM 
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OPTIMUM  PERFORMANCE  OF  A  BUFFERED  SHOCK  TUBE  WITH  A 
URGE  AREA  CHANGE  AT  THE  SECOND  DIAPHRAGM 
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Figure  25  WAVE  PATTERN  IN  A  SIMPLE  SHOCK  TUBE  OF  FINITE  LENGTH 
FOR  THE  PRODUCTION  OF  MAXIMUM  TESTING  TIME 
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Figure  22  WAVE  DIAGRAM  FOR  A  BUFFERED  SHOCK  TUBE  OF  FINITE  LENGTH 
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Figure  23  continued 
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PRESSURE  REQUIRED  IN  DRIVEN  SECTION  TO  MAINTAIN  A  GIVEN  SHOCK 
MACH  NUMBER  AS  DRIVEN  SECTION  DIAMETER  IS  INCREASED 
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Figure  25  WAVE  DIAGRAM  OF  SHOCK  TUBE  WITH  AREA  EXPANSION  AT 
SECOND  DIAPHRAGM 
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Figure  26  WAVE  DIAGRAM  FOR  A  BUFFERED  SHOCK  TUBE  WITH  A 

CONVERGING-DIVERGING  SECTION  AT  THE  SECOND  DIAPHRAGM 
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(jure  27  EFFECT  OF  A  THROAT  AT  THE  SECOND  DIAPHRAGM  ON  THE  PERFORMANCE  OF  A 
BUFFERED  TUBE  WITH  LARGER  DRIVEN  SECTION  THAN  DRIVER 
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EFFECT  OF  VARYING  BUFFER  DIAMETER  ON  THE  PERFORMANCE  OF  A 
BUFFERED  SHOCK  TUBE  WITH  A  LARGER  DRIVEN  SECTION  THAN  DRI 
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Figure  29  HUGONIOT  CURVE  FOR  DETONATIONS  AND  DEFLAGRATIONS 
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